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We construct an evolution equation for the invariant-mass distribution of light-quark and gluon 
jets in the framework of QCD resummation. The solution of the evolution equation exhibits a 
behavior consistent with Tevatron CDF data: the jet distribution vanishes in the small invariant- 
mass limit, and its peak moves toward the high invariant-mass region with the jet energy. We also 
construct an evolution equation for the energy profile of the light-quark and gluon jets in the similar 
framework. The solution shows that the energy accumulates faster within a light-quark jet cone 
than within a gluon jet cone. The jet energy profile convoluted with hard scattering and parton 
distribution functions matches well with the Tevatron CDF and the large-hadron-collider (LHC) 
CMS data. Moreover, comparison with the CDF and CMS data implies that jets with large (small) 
transverse momentum are mainly composed of the light-quark (gluon) jets. At last, we discuss the 
' application of the above solutions for the light-particle jets to the identification of highly-boosted 

, heavy particles produced at LHC. 
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I. INTRODUCTION 

It is known that a top quark produced almost at rest at the Tevatron can be identified by measuring isolated 
jets from its decay. However, this strategy does not work for identifying a highly-boosted top quark produced at 



quark 

the Large Hadron Collider (LHC). It has been observed that an ordinary high-energy QCD jet [l|, can have an 
invariant mass close to the top quark mass. A highly-boosted top quark [3-6], producing only a single jet, is then 
difficult to be distinguished from a QCD jet. This difficulty also appears in the identification of a highly-boosted 
new-physics resonance decaying into standard-model (SM) particles, or Higgs boson decaying into a bottom-quark 
\q ' pair [7|, |8j. Hence, additional information needs to be extracted from jet internal structures in order to improve the 
jet identification at the LHC. The quantity, called planar flow [9(, has been proposed for this purpose, which utilizes 
the geometrical shape of a jet: a QCD jet with large invariant mass mainly involves one-to-two splitting, so it leaves 
a linear energy deposition in a detector. A top-quark jet, proceeding with a weak decay, mainly involves one-to-three 
splitting, so it leaves a planar energy deposition. Measuring this additional information, it has been shown with event 
generators that the top-quark identification can be improved to some extent. Investigations on various observables 
associated with jet substructures using event generators can be found in Refs. 0, llO-pill . For a review on recent 
^ | theoretical progress and the latest experimental results in jet substructures, see Ref. [25|. 

In this paper we shall propose to measure a jet substructure, called the energy profile, which describes the energy 
fraction accumulated in the cone of size r within a jet cone R, with r < R. Its explicit definition is given by [26| 

with the normalization ^(R) = 1, where Pt% is the transverse momentum carried by the particle i in the jet J, and 
ri < r (fj < R) means the flow of the particle i into the jet cone r (R). Different types of jets are expected to exhibit 
different energy profiles. For example, a light-quark jet is narrower than a gluon jet; that is, energy is accumulated 
faster with r in a light-quark jet than in a gluon jet. A heavy-particle jet certainly has a distinct energy profile, 
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which will be studied in a forthcoming paper. The importance of higher-order corrections and their resummation for 
studying a jet energy profile have been first emphasized in (27j . The invariant mass distribution of a single jet has 
also been analyzed in [28| as part of a calculation of threshold effects in dijet cross section. In this work we shall 
apply the perturbative QCD (pQCD) resummation technique (29j . which is extended from the Collins-Soper-Sterman 
resummation formalism [30], to this jet substructure. An alternative approach based on the soft-collinear effective 
theory (SCET) and its application to jet production at an electron-positron collider can be found in Refs. [3ll - l33l |. 

We first derive an evolution equation for the distribution of jet invariant mass Mj, starting with the definitions of 
a light quark jet and of a gluon jet with the four momentum Pj [9|, |3J] . The definition of a jet function contains a 
Wilson line along the light cone, which collects gluons collimated to the light parent particle and emitted from other 
parts of a hadron-hadron scattering process. To perform the resummation, we vary the Wilson line into an arbitrary 
direction tt/' with n 2 ^ (35[. The jet function must depend on Pj and n M through the invariants Pj = Mj and Pj ■ n 
which are related to the jet transverse momentum P T = yj (Pj) 2 — Mj, and n 2 . When Mj approaches zero, the phase 
space of real radiation is strongly constrained, so the associated infrared enhancement does not cancel completely that 
in virtual correction. The infrared enhancement then generates the double logarithms of the ratio (Pj ■ n) 2 / (Mjn 2 ) , 
and the variation of n turns into the variation of Mj. All the different choices of the vector n are equivalent in the 
viewpoint of collecting the collinear divergences associated with the jet. Therefore, the effect from varying n does not 
involve the collinear divergences, which can then be factorized out of the jet, leading to an evolution equation in n 
for the jet function. 

The evolution equation for the jet function is constructed in the Mellin N space, i.e., the space conjugate to 
Mj/(RPt), through which the dependence on the jet cone size R is introduced. Solving the evolution equation, 
we derive the jet function in N as a result of the all-order summation of the double logarithms In 2 N. An inverse 
transformation is then implemented to bring the distribution back to the Mj space. At this step, a nonperturbative 
contribution in the large N region is included to avoid the Landau pole of the running coupling constant and to 
phenomenologically parameterize effects from hadronization and underlying events. This contribution modifies the 
behavior of the jet function at small Mj, but not the behavior at large Mj. It will be shown that our resummation 
results for the jet distribution are consistent with the Tevatron CDF data [36[ . We also observe that a gluon jet has a 
higher invariant mass and a broader distribution due to stronger radiation caused by the larger color factor Ca = 3, 
compared to Cf = 4/3 for a light-quark jet. 

The QCD resummation formula is then extended to the jet energy functions for a light quark jet and for a gluon 
jet, whose definitions are similar to the jet functions. They also contain the Wilson lines along the light cone, which 
collect gluons emitted from other parts of a collision process and collimated to the parent particles. The difference is 
that a step function fcjy©(r — r,-) is associated with each final-state particle i in the smaller jet cone r, where kix and 
are the transverse momentum and the radial distance of the particle i with respect to the jet axis. When r approaches 
zero, the phase space of real radiation is strongly constrained, so the associated infrared enhancement does not cancel 
completely that in virtual correction, which then generates the double logarithms of the ratio (Pj ■ n) 2 /(n 2 r 2 ). The 
derivation of the evolution equation for the jet energy function is basically the same as that for the jet function, and 
the variation of n turns into the variation of r in this case. Because we shall consider the energy profile with the 
jet invariant mass being integrated over, the nonperturbative contribution is not relevant in predicting the jet energy 
profile. The obtained jet energy function allows us to calculate the energy profile 4"(r) in Eq. (JTJ). It will be shown 
that our resummation results for are in agreement with the Tevatron CDF [26[ and LHC CMS [37j data. We 
also observe that a light-quark jet is narrower than a gluon jet, and that jets with high (low) transverse momentum 
are dominated by light-quark (gluon) jets in hadron collisions. 

The above formalism is applicable to the study of a highly boosted heavy particle, with the associated collinear 
radiation being factorized into a heavy-particle jet function. The resultant definition is similar to the light-particle 
jet function, except that the light- particle field is replaced by the heavy-particle field. We then lower the scale to 
the heavy-particle mass tuq, at which jets formed by the light particles, from the heavy-particle decay, are further 
factorized. This step is similar to the conventional heavy-quark expansion, and the factorization of the light-particle 
jet functions holds at leading power of 1/mq. The heavy-particle jet function is thus written as a convolution of a 
heavy-particle kernel, involving specific decay dynamics, and the light-particle jet functions. The former is evaluated 
perturbatively to certain orders of the coupling constant, and results derived in the present work are employed as 
inputs for the latter. Hence, both the heavy-particle jet distribution in invariant mass and the energy profile within a 
heavy-particle jet can be predicted, which will improve the particle identification at LHC. Broad applications of our 
framework to jet physics are expected. 

In Sec. II, we construct the evolution equations for the light-quark and gluon jet functions, and solve them in the 
Mellin space. The treatment of soft gluon contributions to the evolution equations is explained. A nonperturbative 
contribution is introduced into the resummation formula to mimic PYTHIA8.145 (38j predictions in the region of 
small jet invariant mass. After fixing the nonperturbative piece at a given Pt value, the behavior of the jet functions 
in the whole range of invariant mass is derived via the inverse Mellin transformation numerically in Sec. III. It will 
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be shown that our resummation predictions for the jet mass distribution agree well with the CDF data. The same 
formula is extended to calculating the energy profiles of the light-quark and gluon jets in Sec. IV by constructing 
and solving the evolution equations for the jet energy functions. Our resummation predictions are consistent with 
the CDF and CMS data. With the important logarithms being collected, the initial conditions of the jet functions 
and the jet energy functions can be evaluated up to a fixed order. Their ncxt-to-leading order (NLO) expressions are 
presented in Appendices A and C, respectively. The contour choice for the inverse Mellin transformation is discussed 
in Appendix B. Before concluding this section, we note that the non-global logarithms and the clustering effects should 
be also considered, when comparing experimental data and theoretical predictions for the jet mass distribution at the 
next-to- leading- logarithmic (NLL) level, as discussed in Refs. [39l - l4]| . 

II. RESUMMATION FOR JET FUNCTIONS 

In this section we derive the evolution equation for the light-quark and gluon jet functions defined in 34]: 

(27T)" 



J q (Mj,P T ,^,R,fi 2 ) = 



o Avow E Tr {^ lg(°)^ ?)t (^' )l^)^l^n ?) (°o,0)9(0)10)} 



xS(M 2 - M 2 j(Nj, R))5W(e - e(Nj))6(P°j - u(Nj)), 

xS(M 2 j - Mj(Nj,R))S^(e - e(Nj))S(P° T - lo(Nj)), (2) 

where \Nj) denotes the final state with Nj particles within the cone of size R centered in the direction of the unit vector 
e = (0, 1,0,0), Mj(Nj, R) (ui(Nj)) is the invariant mass (total energy) of all Nj particles, and /i is the factorization 
scale. The above jet functions absorb the collinear divergences from all-order radiative corrections associated with 
the energetic light jet of momentum Pj = PjV^, where Pj is the jet energy, and = 0,0) is a 4-vector with 
P = y/l — (Mj/Pj) 2 . The coefficients in Eq. ([5]) have been chosen such that the lowest-order (LO) jet functions 
are equal to S(Mj) in perturbative expansion. The definition of the jet function in Eq. ([2]) contains a Wilson line, 
which collects gluons radiated from either initial states or other final states of a hadron-hadron scattering process, and 
collimated to the light-quark (or gluon) jet. Gluon exchanges between the quark fields q (or the gluon fields F' 71 ' and 
F£) correspond to final-state radiation. Both initial-state and final-state radiations are leading-power effects in the 
factorization theorem, and have been included in the jet function definition. However, the contribution from multiple 
parton interaction, which is regarded as being higher-power, is not included. Nevertheless, it still makes sense to 
compare predictions for jet observables based on Eq. ([2]) at the current leading-power accuracy with experimental 
data . 

The Wilson line represents the path-ordered exponential 



W„(oo,Q) = Pexp 



ig s (fi 2 ) / dzn ■ A(zn) 



o 



(3) 



where the gauge field denotes A = A a t a with t a being the gauge group generators in the fundamental (adjoint) 
representation for the light-quark (gluon) jet function, and g s ([i 2 ) is the QCD strong coupling at the energy scale /i. 
As explained in the Introduction, the original Wilson line vector £ = (1,-1,0,0) [34j can be replaced by the arbitrary 
vector n, while the spin projector in the light-quark jet, cf. Eq.([5]), remains unchanged. The scale invariance of 
Eq. ([3]) in n guarantees that the jet function depends on the ratio 

= 4(vn) 2 

~ R 2 \n 2 \ ' 1 ' 

where the dependence on R is inspired by the logarithms observed in the NLO jet function. We then vary n by 
considering the derivative [35[ of the jet function Jf : 

n 2 d 



v a — J f {M 2 ,P T ,v 2 1 R^ 2 ) 1 (5) 
V ■ n dn a 

with / = q or g. The n dependence appears only in the Feynman rules for the Wilson line, whose differentiation with 
respect to n a leads to 



-Vn 



d n„ n 2 ( v ■ I \ 1 n,. 

- p (6) 



v ■ n dn a n ■ I v • n \n ■ I I n ■ I n ■ I 
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The special vertex defined in the above expression suppresses the collinear region of the loop momentum I that 
flows through the special vertex: if I is parallel to Pj, i.e., to v, the contribution from the first term is down by the 
ratio Mj/Pp. The second term u M also gives a power-suppressed contribution, after being contracted with a vertex 
in Jf, in which all momenta are mainly parallel to Pj, Hence, the leading regions of I are soft and ultraviolet, but 
not collinear. 





FIG. 2: Factorization of the LO real soft kernel. 



To obtain the leading logarithms (LL), the special vertex must appear at the outermost end of the Wilson line 
(nearest the final-state cut) as shown in Fig. HJa). If the special vertex does not appear at the outermost end, the 
gluons emitted after the differentiated gluon must be soft too. Otherwise, their finite momenta will regularize the 
soft divergence associated with the differentiated gluon. In this case we will have more soft gluons, namely, a soft 
divergence at higher orders in the coupling constant, which corresponds to a subleading logarithm. To collect the 
LL in Fig. [TJ the replacement g**" — > Pjl v / (Pj ■ I) [42| is employed for the metric tensor of the differentiated gluon, 
where the vertex with the Lorentz index fi is located on the Wilson line, and the vertex v on a line in the jet function. 
We explain this replacement by assuming that Pj is in the plus direction for convenience. Then the component 
among gV leads to the leading contribution. The + superscript is represented by the largest component Pj~ of Pj 
in the replacement. The components V are arbitrary, but only l~ is selected when V is contracted with a vertex in 
the jet function, which is dominated by the momentum flow along Pj. Applying the Ward identity to the sum over 
all possible attachments of V (42[, we factorize the differentiated gluon into the virtual soft kernel as displayed 
in Fig. [TJ The factorization of the real soft kernel at LO is depicted in Fig. [2] The LO soft kernel is then 
written as the sum of the above two diagrams, i.e., 

To produce a LO ultraviolet divergence, the special vertex must appear at the innermost end of the Wilson line, 
and the differentiated gluon forms a loop correction to the quark- Wilson-line vertex as shown in Fig. [3J If this is not 
the case, we will have more off-shell lines, namely, a higher-order ultraviolet divergence, which leads to a subleading 
logarithm. The LO differentiated gluon can be factorized trivially by performing the Fierz transformation of the 
fermion flow, 

hjhk = + ^(75) l fc(7s)ii + ^(l a )zkh a )ij + ^hsla^k^^i] + ^{o a p)ik{(J a? )ij, (7) 
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with / being the identity matrix, and a a p = i[y a , 7^]/2. The first and last terms contribute in the combined structure 



lull 



ijJ-lk 



(8) 



where the vector £ lies on the light cone and satisfies £ • £ = 1. The identity matrix In- in Eq. © goes into the trace 
for the jet function. The matrix (£ / £);j/4 then leads to the loop integral for the hard kernel G^ in Fig. [3] 

The jet transverse momentum, the jet invariant mass, and the jet cone, under the factorization of the virtual 
differentiated gluons, remain as P F , Mj and R, respectively. The jet momentum and the jet cone arc not modified 
by the soft real correction, but the jet invariant mass squared Mj, regarded as a small scale, is modified into 
(Pj — I) 2 — M 2 — 2Pj ■ I. For the light-quark jet function, we then arrive at the differential equation 



-J q (M 2 ,P T ,v 2 ,R,^ 2 ) ^2(G + K) ® J q (M 2 ,P T ,y 2 ,R^ 2 ), 



(9) 



where the hard correction, the virtual soft correction, and the real soft correction to the NLO evolution kernels are 
written as 



Ki 1 



® J„ 



ig s C F ^ 
-ig 2 C F ^ £ 



'I 



(2tt) 4 - £ (n ■ I + ie)(l 2 + ie) 
d*- e J n ■ P„ 



4 



P 



Jv 



(27r) 4 - e [n-l + ie){Pj ■ I - ie)(l 2 + ie) 



(Pj - I) 2 + IE 

SK, 



P.T-1 



6G, 



g 2 s Cj 



(Pi 



n-Pj 



(2ir) 4 (n ■ I + ie)(Pj ■ I - ie) 



—2nS(l 2 )J q (Mj - 2Pj ■ I, P T , v 2 ,R, n 2 ), 



(10) 

(11) 
(12) 



respectively. The first term in the parentheses of Eq. (fT0| is free of ultraviolet divergence, and the second term, repre- 
senting the soft subtraction —K^ to avoid double counting of the soft contribution, contains ultraviolet divergence. 
As adding G^> and K^ together, their ultraviolet divergences cancel. K^ in Eq. (fT2j) is ultraviolet finite, so the 
kernel G + K = G + K v + K r is independent of renormalization scale fjf. In our regularization scheme, the additive 
counterterms SG and SK are chosen as 



SG = —C F 
2tt 

= -SK, 



- +ln(47rC 2 V) - lE 



(13) 



where a s = g 2 /4ir, je is the Euler constant, and the arbitrary constant C2 can be varied to estimate subleading 
logarithmic corrections to our formula. 

The trace in Eq. ([TU|) indicates that the v v term in the special vertex n v gives a contribution suppressed by Mj/P^, 
as compared to the contribution from the n v term. Equation IjlOp then reduces to 



G« = ig 2 s C Ff i" 



Z7T 



Pj-n 



(2*y 



fi,(f>j- fiPj-l Pj-n 



(n-l) 2 (Pj ~l) 2 l 2 (n-l) 2 l 2 



6G, 



: (c 2 u 2 rp t ) 2 

m - 1 



(14) 
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The virtual soft correction in Eq. (ITT1) gives 



v-(l) -2^, It 2 I d 4 e l 2Pj ■ I 



CF^^TTT-a^ (15) 



in which the infrared regulator A 2 will be taken to be zero eventually. 

It is more convenient to perform the resummation in the conjugate space via the Mellin transformation. The reason 
becomes evident as comparing the convolutions of the virtual and real soft corrections with the LO jet function: the 
former leads to <g> J (0) = Ki^8(Mj), while the latter leads to 



= < 16 » 



If transforming the above results into the Mellin space, the infrared divergences from Mj — > in the virtual and real 
soft corrections cancel explicitly. Therefore, we introduce the Mellin transformation 



l 

N-l r („ D _ 2 D ,.2\ 



J q {N,P T ,v\R,^) = dx(l -x) iv ^ i J 9 (a;,P T ,^,i?,^), (17) 
Jo 

x = M 2 j/(RPt) 2 being the dimensionless variable. The convolution in Eq. (| 12[) is converted into a product 

/ dxil-x) 1 *' 1 ^ ® J q = KP(N)J q (N,P T ,u 2 ,R,n 2 ), (18) 
Jo 



with the definition 



mm = [ mi - .)»-' / ^ (n . l+ % J Xu*) i(l2)i (• - 2 A (1 - cos9) ) ■ (19) 

To derive the above expression, we have made the small-mass approximation 1 — ft cos 9 ~ 1 — cos#, and inserted the 
identities J dz8(z — 2|1|(1 — cos 6) J (RPx)) — 1 and J dy5(x — y — z) = 1. The approximation 1 — x = 1 — y — z « 
(1 — y)(l — z) has been also adopted, which holds in the dominant region with small y and z. 
We compute Eq. (fT9f by splitting it into two pieces 



K?\N) = ^C F j\[(l-z) If - 1 -l] J (2 t)3 (n . "+ ic)^ * (* ~ 2 R% (1 ~ COSg) ) Q(E ~ g) 

+^ / , ,S;2 nfl , , ia ^ a a )gf g -2j|-(i- C oBg)V (20) 



II 



(27r) 3 (n-/ + ie) 2 (2Pj-/ + A 2 ) v ; V R p T 



where the infrared regulator A 2 has been neglected in the first term, because of the absence of the infrared divergence 
from z — > 0. Since the gluon momentum is finite in the first term, we require that its angle can not exceed the cone size 
R by including the step function Q(R — 8), which then brings the R dependence into our resummation formula. The 
soft effect dominates in the second term, so there is no need to constrain the range of the angle 9. A straightforward 
calculation leads to 

KW(N) = ^C F lnE§, (21) 
7r TV A -2 



with TV = A r exp(7e). Combining Eqs. (|15|) and (|2~T|) . we obtain 

7T 



CiRPt , . C*2 
In — = h In — 

Nfi' Ci 



(22) 



where = and the dependence on the infrared regulator A 2 has disappeared. Furthermore, an arbitrary 

constant C\ has been introduced to estimate subleading logarithmic corrections to our formula. 
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Solving the renormalization-group (RG) equations, 



with the cusp anomalous dimension 



(23) 



(24) 



we derive 



A" 



— ,<W 



K (l,a a 
Cn 



G 



C 2 v 2 RPt 



/2- 



C 2 ii 2 P 2 



N 2 

C 2 R 2 P 2 
N 2 



In 



+ G(l,a s (C 2 Vi? 2 P|)) - / 

./c 

a., (C 2 Vi? 2 P|) 



C 2 u RPt 



dp,' 



dRP T /N M 



-Ajc(a.(/x")), 



^ + fr^- 1 »-''>- 

Cl 27T 



Ci/iV 



do; 



\ K {a s {u 2 R z P^)) 



(25) 



With the large logarithms being removed, the LO expression for the initial condition K(l,a s ) + G(l,a s ) of the RG 
evolution has been inserted into the last line. The cusp anomalous dimension Aa' is process independent, and given, 
up to two loops, by 



7T F ' 2 



r I 67 ^ 



— nt 
9 ' 



(26) 



for a light quark jet, where rif denotes the number of active light-quark flavors. 

After organizing the large logarithms in the kernels, we solve the differential equation 

tl 2 d — d — 
v a — j q {N,P T ,v 2 ,R^ 2 ) = 2v 2 — J q {N,P T ,v 2 ,R,Li 2 ) 



\ Nfj,' ' ^ ' 



C 2 v 2 RPt , /2 . 
G\ ; ! a s\p- , 



J q (N,P T ,v 2 ,R,fi 2 ). 



(27) 



The strategy is to evolve v 2 from the low value z-f n = Cx/iC^N) to the large value v\ = 1, corresponding to the specific 
choices n = n in = (1, (4C 2 A - C*iA 2 )/(4C* 2 A + C*iA 2 ),0,0) and n = n R = (1, (4 - A 2 )/(4 + R 2 ),0,0), respectively. 
The former defines the initial condition of the jet function, which can be evaluated at a given fixed order, because of 
the vanishing of the logarithm ln(C 2 ^ 2 -/V/Ci). The latter defines the all-order jet function with the large logarithms 
being factorized and organized. Since the jet function collects the soft and collinear radiations, which mainly occur 
at a lower scale, /j 2 should take a value of 0(R 2 P^/ N). This choice introduces an additional single logarithm, that 



needs to be summed to all orders by a RG evolution equation in \i. To achieve it, we set \j? 
will be elaborated in Appendix A. The solution to Eq. ([2"7)l is derived as 

J q {N,P T ,vl,R) = J q (N,P T ,v 2 n ,R)eMS q (N,P T ,R)}, 

with the Sudakov exponent 

r c * dy 



0{R z Pj,/{Nv 2 )), which 



(28) 



S q (N,P T ,R) 



Ci/jv y yjd/N 



du 



X K {a s {u 2 R z P 2 )) 



2tt 



a s (y 2 R 2 P 2 ) 



C F n ( C{R 2 P 2 
N 2 



In 



9i 

Ci 



(29) 



ft is noted that the 7? dependence appears in the single logarithmic term of the Sudakov exponent. 
We further evolve a s from the scale C\RPt/N to yRPr in the last term of Eq. (|29|) . 



C F 



fC 2 R 2 P 2 \ 






V N 2 J 


7T 


/ 



r a B (C 2 R 2 P 2 /N 2 ) 



C, 



yRPr 



dfl 



dRP T /N M 



2/3(a s ( M 2 )) 



da s +a s (y 2 R 2 P%) 
a s {y 2 R 2 P 2 ) 



(30) 
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and expand the QCD Beta function up to 0{a 2 ), f3 = -(P /A){a s /n;) 2 with (3 = 11 - 2n//3 [43]. Inserting Eq. <J30]) 
into Eq. (|29p . and applying the integration by part, the exponent is rewritten as 



C2 rfy 



S q (N,P T ,R) = 
with the anomalous dimensions 

A q = C F ^ + \C F (^) 



A 9 K(y 2 i? 2 P|))ln(^ ) +^ / (n,i,;-/?-7';)) <► . 



1 18 6/ 9 ' M0 C*i 



(31) 



The Sudakov exponent for the gluon jet function can be derived in a similar way 

r ' 1 dy ( . , , o„o„oss , ( Q 

ICi/N V 



S g (N,P T ,R) 



^-{A g {a s {y 2 R 2 P 2 ))hx[^ 



B g (a s (y 2 R 2 P 2 )) 



(32) 



(33) 



where the anomalous dimension A g (B g ) is obtained by substituting Ca for Cf in A q (B q ). In this work the NLL 
terms have been included into the resummation by adopting Af at two- loop level and Bf at one-loop level. Although 
the numerical evaluation of the Sudakov integral induces some next-to-next-to-leading logarithmic (NNLL) terms, 
the inclusion of the complete NNLL terms demands higher-order contributions to Af and Bf. Hence, we shall refer 
our resummation formalism presented here as one with the NLL accuracy. Finally, it is noted that the non-global 
logarithms discussed in Refs. [39l - l4ll | are not included in our resummation formalism for the jet function definition in 
Eq. ©. 

We evaluate the initial conditions of the Sudakov evolution for the light-quark and gluon jet functions up to NLO 
in Appendix A, and confirm that the large logarithms In TV do not appear in these initial conditions as v 2 — v 2 n ; 
namely, they have been collected into the Sudakov exponents. We note that the quark-loop contribution to the gluon 
jet function, which carries a different color factor, has to be handled separately as shown in the next section. The 
resummation formulas for the light-quark and gluon jets are summarized, in the Mellin space, as 



J q {N,P T ,R) 
J g (N,Fr,R) 



1 



R 2 P 2 



1 + ^L as (C 2 R 2 P 2 ) 



1, d 1, 2 Cx 

2 ln c^-2 ln c; 



xS q (N,P T ,R), 

R 2 P 2 { n 
xS g (N,P T ,R), 



1 + ^A as (c 2 R 2 P 2 



1, C t 



In 



c 2 



1 Cg£ 

4 C 2 



lln^ 
12 C 2 



1 



9 



2 lE 4 8 



(34) 
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IE 



TV" 

T 



^ ln2 - 3 ) + 3^ 



(35) 



Here the third arbitrary constant C3 has been introduced through the choice of the renormalization scale [i for the 
initial conditions, which denotes another source of theoretical uncertainty in our formalism. 



III. NUMERICAL ANALYSIS FOR JET FUNCTIONS 



In this section we compare our predictions for jet mass distribution to the experimental data from the Tevatron and 
the LHC. As x — Mj / (RPt) 2 — > 0, all moments in N are equally weighted, since the suppression factor (1 — .t)^ -1 
is not effective. The terms containing In AT, being the dominant ones, have been summed to all orders in a s , so the 
predictions from Eqs. (|34"|) and (|35|) are supposed to be reliable at small x. However, the running coupling constant a s , 
evaluated at the soft scale RPt/N, increases with N, and the expansion parameter a s laN may become much larger 
than order unity. In this region a perturbative calculation is not adequate and contributions from nonperturbative 
physics need to be included. Furthermore, the complex argument \i — yRPr of a s {^ 2 ) in Eqs. (|3"Tt and (f3"5]) tends 
to be small in magnitude at large N, even lower than the Landau pole scale. Therefore, in our numerical analysis 
we introduce a critical scale /i c to avoid the Landau pole, below which the running coupling is frozen to the constant 
value a s (fj, 2 ). For an explicit treatment of a s (fi 2 ), see Appendix B. As x grows gradually, the large- iV moments are 
suppressed by (1 — x)^ -1 , and the resummation effects together with the nonperturbative inputs become less crucial. 
A fixed-order evaluation is then more reliable at large x, where Eqs. (j34f and (|35"T) are expected to coincide with the 
NLO jet mass distributions, cf. Appendix A. 
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In this work the following nonperturbative correction is implemented into the Sudakov exponent in the N space 



Sf p (N,P T ,R) 



N 2 Ql 
R 2 P 2 



(C/ao IniV + ai) + C/a2 



RP T ' 



(36) 



with Q = 1 GeV and Cf = Cf(Ca) for the light-quark (gluon) jet function. The first two terms proportional 
to N 2 Q^/P T are similar to the singular terms in the nonperturbative contributions to the transverse- momentum 
resummation [30l l44l |4o| and threshold resummation (46j formalisms. The last term, being a power correction [47j j . 
can be obtained from the asymptotic behavior of the Sudakov exponent. The powers in NQq/Pt indicate that 
the nonperturbative effects are significant only in the extremely large N region. We determine the nonperturbative 
parameters cto, ct\ and «2 from fits to PYTHIA8.145 [38[ predictions associated with SpartyJet [H[ for the light-quark 
and gluon jets, separately. The resummation formulas including the nonperturbative inputs are then written as 



Jf ES (N,P T ,R) = J q (N,P T ,R)exp[S^ p (N,P T ,R)}, 

JJ es (N,Pt,R) = J g (N,P T ,R)eMSr(N,P T ,R)] + 3^0. (l ~ * G 



(37) 
(38) 



where the quark- loop contribution proportional to the flavor number n/ has been added as the second term on the 
right-hand side of Eq. (|38|) . Note that this contribution does not contain the large logarithm IniV as /j, 2 ~ 0(R 2 P^/N), 
at which the final conditions of the jet functions are defined, so it is not organized into the resummation formula. 
The inverse Mellin transformation of the above expressions leads to 



jf* s (M},P T ,R) 



1 



dN(l 



- N jf ES (N,P 7 



R). 



(39) 



An appropriate contour C extending to infinity in the complex N plane needs to be chosen for the numerical inverse 
transformation, which is specified in Appendix B. 

As stated before, hard radiation is important at large Mj, although the probability of having a jet with large mass 
decreases quickly as Mj increases. To describe the distribution at large Mj, we further perform the matching between 
the resummation and NLO results via 



jNLL/NLO^ R) = JRES (M 2 ; p y _ R) + [jW^Mj, P t , R) - J^ R ' asym (M 2 j , P T , R) 

J g NLL/NLO (Mj, P T , R) = Jf E3 (Mj, P T , R) + \j { g 1)R (M], P T , R) - J^ R ' asym (Mj, P T , R) 



(40) 



where Jt is the contribution from the NLO real emissions, jW fl,asym denotes its asymptotic expression in the 
Mj —> limit, i.e., the so-called "singular piece" [2{J. The inclusion of the "regular piece", i.e., the term in the square 
brackets on the right-hand side of Eq. (|4T)|) . warrants that the expansion of jJ LL / NLO U p ^ NLO coincides with the 
complete NLO QCD predictions of the jet functions. We note that the regular piece of the quark-loop contribution 



to the gluon jet function has been included into J, 



(i)R 



J, 



(l)-R.asym 



, cf. Appendix A. 
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FIG. 4: Quark (left) and gluon (right) jet mass distributions with S NF (solid lines) and without 5* NP (dotted lines) for Pt = 600 
GeV and R = 0.7. 



To be compared with the normalized jet mass distribution, we convolute Eq. (|40[) with the parton-level differential 
cross section daf/dPr evaluated at the renormalization scale jj, = C^RPt, the same as the initial scale in Eqs. (|3"4")l 
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and (|35p . yielding the factorization formula 



1 da 



1 



E / dP T ^(Mj,P T )J^°(Mj,P T ,R), 



(41) 



a dM 2 T 



a 



where a = j(da/dMj)dM 2 is the integrated jet cross section. We adopt the default choice C\ = exp(7s), C2 = 
exp(— 7b), C3 = 1, and /i c = 0.3 GeV, and include the nonperturbative contributions in fits to PYTHIA predictions 
for the jet distributions with Pt — 600 GeV and R — 0.7. It is found that the nonperturbative parameter set 
ao = —0.35, ai — 0.50 (a\ = —4.59), and a2 = —1.66 leads to a reasonably good fit to the light-quark (gluon) jet. It 
is also observed that the quark-loop contribution to the gluon jet function is negligible. 

The quark and gluon jet mass distributions depicted in Fig. [4] indicate that including S" NP shifts their peak positions 
toward the larger jet mass region, and suppresses (enhances) the peak height of the quark (gluon) jet distribution. As 
stated in the Introduction, the nonperturbative contribution does not modify the behavior of the jet functions at large 
Mj. Given the nonperturbative parameters, we predict the jet mass distributions at any arbitrary value of collider 
energy jet energy Pt and jet cone size R. The resummation predictions for the normalized light-quark and gluon 
jet mass distributions as functions of Mj /(RPt) for R = 0.4, 0.5, 0.6 and 0.7 with RPt = 280 GeV are presented 
in Fig. [5J It has been found in (4(| that the NLO jet mass is remarkably well described by the simple rule-of-thumb 
Mj ~ 0.2RPt- However, Fig. [5] shows that not only the average jet mass but also the shapes of the light-quark 
and gluon jet mass distributions almost remain the same, when we vary the jet cone R with RPt being fixed. This 
behavior is attributed to the fact that each component of the resummation formula, including the Sudakov factors in 
Eqs. (|3"Tj) and (|33[) . the initial conditions in Eqs. (fM]l and (l3"5j) . and the nonperturbative contributions in Eq. (|36l) . 
depends only on the scale RPt- The scaling behavior is violated when the jet mass is large enough (Mj /(RPt) > 0.7), 
as indicated in Fig. [5j Nevertheless, the probability to find a jet with such a large mass is low. We also note that 
the jet mass distribution as a function of Mj/(RPt) is relatively independent of the collider energy y/~S, except 
that for substantially larger momenta the reduced phase space will lead to smaller predicted jet masses at the same 
momentum. Furthermore, our formalism also suggests that this conclusion holds for a similar jet (with the same Pt 
and R) produced in any kind of hard scattering processes, such as the associated production of jets with gauge boson 
or Higgs boson. 

Following Eq. (14"TT) . we convolute the light-quark and gluon jet functions with the constituent cross sections of LO 
partonic dijet processes at the Tevatron and the parton distribution functions (PDF) CTEQ6L [50]. Here we have 
neglected the soft gluon contribution [5l|, equivalent to the soft function introduced in the Soft Collinear Effective 



Theory (SCET) [521 ] . which couples the light-particle jet and the partonic processes. The resummation predictions 
for the jet mass distributions at R — 0.4 and R = 0.7 are compared to the Tevatron CDF data [36| in Fig. [5] with 
the kinematic cuts Pt > 400 GeV and the rapidity interval 0.1 < \Y\ < 0.7 . The above data were obtained using 
the midpoint jet algorithm [53j . and the data from the anti-fc t algorithm [54] do not vary much as shown in [36]. The 
consistency of the resummation results with the CDF data is excellent at intermediate Mj. The resummation formula 
describes the shapes and the peak heights of the jet distributions in the small Mj region, but with the peak positions 
being slightly lower than the CDF data. As indicated in (36|, the PDF uncertainties could induce large variation in 
shapes of jet mass distributions around peak positions. The difference from the data in Fig. [5] is within the above 
variation. This is the first time that the pQCD factorization theorem explains the observed jet mass distributions 
successfully. Note that the jet mass distribution, which corresponds to the angularity distribution with a — [3lj . 
cannot be well described in the SCET formalism. In Fig. [7] we display the resummation predictions for the jet mass 
distributions at the Tevatron with R = 0.3 and at the LHC with R = 0.7, which can be tested by Tevatron data and 
LHC experiments. 



We define the jet energy functions Jf (Mj, Pt^v 2 , R,r) with / = q(g) denoting the light-quark (gluon), which 

describe the energy accumulation within the cone of size r < R. The definition is chosen, such that jf — Pr^(Mj) 
at LO. In this section we will study the energy profile of a light-particle jet in the framework of QCD resummation at 
leading power of r. The Feynman rules for Jj are similar to those for the jet functions Jf at each order of a s , except 
that a sum of the step functions J2i k®Q(r — 9{) is inserted, where (6i) is the energy (the angle with respect to the 
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FIG. 5: Resummation results for the light-quark (upper) and gluon (lower) jet mass distributions as functions of Mj /(RPt) 
including the nonperturbative contributions for R = 0.4, 0.5, 0.6 and 0.7 with RPt = 280 GeV. The ratios relative to the 
predictions for R — 0.7 are also shown. 
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FIG. 6: Comparison of resummation predictions for the jet mass distribution to Tevatron CDF data with the kinematic cuts 
Pt > 400 GeV and 0.1 < \Y\ < 0.7 at R — 0.4 and R = 0.7. The inset shows the detailed comparison in large jet mass region. 
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FIG. 7: Resummation predictions for the jet mass distribution for Tevatron and LHC. The inset shows the detailed behaviors 
in large jet mass region. 



jet axis) of the final-state particle i. For example, the jet energy functions Jj are expressed, at NLO, as 

Mm § lwhwk°^ (r - e - )+ k ° e(T - 

xTr[mq(0)W^H^,0)\p,a-k 7 \)(k,X;p,a\W^H^,0)q(0)\0)} 
xS(Mj - (p + k) 2 )6^(e - e p+k )S(P°j -p°- k°), 

x (01^^(0)^^(00, 0)\p, o- k, A) (k, X;p, a\W^ (oo, 0)J*(0)&|0) 
xS(Mj - (p + kf)8^\e - e p+k )S(P°j - p° - k°), (42) 

where the expansion of the Wilson links in a s is understood. As shown in the previous section, the quark-loop 
contribution to the gluon jet function is not important, cf. Eq. (|38j) . with a proper choice of the factorization scale /i 
in the resummation calculation. Hence, the quark-loop contribution to the energy profile of the gluon jet can also be 
ignored with an appropriate choice of fx. 

When r approaches zero, the phase space of real radiation is strongly constrained, so infrared enhancement does not 
cancel exactly with that in virtual contribution and results in large logarithms, e.g., a s In 2 r. An evolution equation 
for summing these logarithms to all orders in a s in the jet energy functions can be constructed, whose derivation 
is similar to that for the jet functions discussed in Sec. II: the variation of the Wilson line direction introduces the 
same special vertex in the differentiated jet energy functions. The virtual gluons emitted from the special vertex are 
factorized into the same hard kernel and the same virtual soft kernel Kv . For example, their expressions for 
the light-quark jet energy function are given by Eqs. (fl4)) and (fl5)) . respectively. For the real soft gluon emitted 
from the special vertex, we split the sum of the step functions into 

J2 k°Q(r - Bi) = £ k°G(r B v ) + l°Q(r - 9), (43) 

i %' 

in which J\, means a summation over final-state particles with the real soft gluon being excluded. The first term in 
Eq. dUl) gives 

*™»Jf = fa J £L . ; - * K g — ^8 (, - flfeg) JfM - aft . ,, ft, A R, r). (44) 

Because of the real soft gluon emission with the polar angle 9, the jet axis of the rest of particles, described by 
on the right-hand side of the above expression, inclines by an angle |l|sin#/Py with respect to the jet momentum 



J^(Mj,P T ,u 2 ,R,r,fj, 2 ) = 



J^(M 2 j,P T ,v 2 ,R,r,n 2 ) = 
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Pj . The step function in Eq. (l44l) imposes a phase-space constraint on the real soft gluon emission, such that the jet 
axis of the rest of particles cannot move outside of the jet cone r. Applying the Mellin transformation with respect 
to x = Mj/(RP T ) 2 , we have 

/ dxil-x)"- 1 ^® Jf = KW(N)jE(N,P T ,v 2 ,R,r), (45) 
Jo 

with the definition 

Kr {N) = 9 s Cf / dz{l-z) N 



(2tt) 3 (n-l + ie) 2 (2Pj ■ I + X 2 ) 



x4( ^(,_A (1 _ Mlfl )e(r-lte»). (46, 
The second term in Eq. (1431) leads to 

whose Mellin transformation gives 

.1 

4b(1 - a?)* -1 !^ <8 J q = KP(N)J q {N,P T y ,R), (48) 



with the definition 



*<-><*, . *c, £ *(. - / ^f£>* (.-*(.- cos,)) e (f - ,„ 



(49) 



Strictly speaking, the energy |1| of a real gluon cannot approach infinity, so the step function at the end of the above 
expression has been introduced. Working out the above integration, we obtain 

RP(N) = ^C F ±- [ dcosO— — —r . RPt „. (1 - cos^ 6) 9(r - 0), (50) 

e K > 2tt N J (n° - n x cos8) 2 2(1 - cosfl) v ' v h y 1 

which is down by 1 — cos^ r and negligible in the small r region. This result is attributed to the suppression of the 
second term in Eq. (|43|) by soft Hence, this piece will not be considered from now on. 

The jet energy profiles are measured by summing over all jet invariant masses in experiments. Therefore, we 
perform a corresponding analysis with the M 2 dependence being integrated out of the jet energy profiles, namely, by 
considering only the N = 1 moment. A straightforward computation leads Eq. (J^HJ) to 

K?\i) = |^m^^!, (5i) 

where the infrared regulator A 2 will be taken to be zero eventually, and v 2 is defined as in Eq. Q. Using the same 
countcrterm, Eqs. (|15l) and (I51|) are combined to form 



= ^ 1) +^ 1) (l) = g^ln^^ + ^C F lng. (52) 



which contains the large single logarithm lnr. 
Solving the RG equation for the kernels, 



M'/yG = X K = -l^'-ffK, (53) 



we derive 



K 



dji 1 ' dfi' 

yPTrC x , 2 \ ^f v 2 c 2 RP T 2 

; — + G } ) 



K (l^s^P^Cl)) +G(l,a s (^C 2 2 R 2 P 2 )) -\ ( C2UPtR ^l XK (a s ^' 2 )), 



^ln^a s {y 2 P 2 r 2 C 2 ) + ^a s (^C 2 2 R 2 P 2 ) - \ f^' ^\ K (a s (u>P 2 )). (54, 
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The light-quark jet energy function then obeys a differential equation similar to Eq. ([9]) : 

,2 



—v a J-jE(l,P T , V 2 ,R, r, ^ 2 ) = 2v 2 -^ Jf (1, P T , v\R, r, 



, , vPrrd fif2 ,\ , ^( v 2 C 2 RP T 2 
K\ } — ,a s (/z ) )+G[ - ,a s (// ) 



Jf (l,iV,^,i2,r,/i 2 ). 



(55) 



A similar equation also holds for describing the energy profile of the gluon jet. As solving these equations, we choose 
the factorization scale /j 2 ~ 0(r 2 P^/(R 2 v 2 )) : so that the quark-loop contribution to the gluon jet energy profile can 
be ignored, for the quark-loop contribution does not contain the large logarithm ln(i? 2 /r 2 ) with this choice of the 
scale. 

The strategy to solve the above equation is to evolve v 2 from the low value v 2 = v 2 n = C 2 r 2 / (C 2 R 2 ) to the large 
value v 2 = v\ = 1, which correspond to the specific choices n — n ln = (1,(4(71 — r 2 C 2 )/(4Cf + r 2 C 2 ), 0,0) and 
n = ne = (1, (4 — i? 2 )/(4 + R 2 ), 0, 0), respectively. The solution of the above equation is given by 



JZ(l,P T ,v 2 ,R, r ) = jZ(l,p T ,vl,R,r)exp[S q (R,r)}, 



(56) 



with the Sudakov exponent 

rC 



S q {R,r) 



dy 
C dy 



dp , C 2 



2tt C 2 



hi ^o. (yP^Cl) + ^a s (,rcWrj.) 



2„„2 id 



Cfyr 



-X K (a s (ioP 2 )) 



2 

^a s (y"C 2 R 2 P 2 ) Q + lng) - i jT ^A g (a>C 2 i? 2 P 2 )) 



(57) 



where the constant C will be fixed below. The Sudakov exponent S g (R, r) for the gluon jet is obtained by substituting 
the color factor Ca for Cp in the above expression. The resummation formulas are summarized as 



Jf(l,P T ,4,R,r) = J?(l,P T ,vl,R,r)e X p[S f {R,r)} 



(58) 



with the subscript / = q or g. The 0(1) constants are chosen as C\ = C2 = 1 and C = exp(5/2) (C = exp(17/6)) in 
order to reproduce the single logarithm a s lnr in the NLO light-quark (gluon) jet energy function. The initial condi- 
tions Jf : (1, Pt ,f 2 n , R,r) of the Sudakov evolution, in the absence of the large logarithms and with the factorization 
scale ix ~ O(Pt), are calculated up to NLO in Appendix D. 
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FIG. 8: Resummation predictions for the energy profiles of the light-quark (solid curve) and gluon (dotted curve) jets with 
= 7 TeV and 80 GeV < Pt < 100 GeV. 

Inserting the solutions in Eq. (1581) into Eq. ([T]) , the jet energy profile is written, in terms of the convolution with 
the parton-level differential cross section, as 



*(r) 



(59) 
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which respects the normalization ^(R) = 1, and vanishes as r — > 0. Note that a jet energy profile, with N = 1, 
is not sensitive to the nonperturbative contribution, so our predictions are free of the nonperturbative parameter 
dependence, in contrast to the case of describing the jet invariant mass distribution, cf. Sec. II. We find that the 
light-quark jet has a narrower energy profile than the gluon jet, as exhibited in Fig. |U for y/S = 7 TeV and the interval 
80 GeV < Pt < 100 GeV of the jet transverse momentum. The broader distribution of the gluon jet results from 
stronger radiations caused by the larger color factor C 'a = 3, compared to Cf = 4/3 for a light-quark jet. 
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FIG. 9: Comparison of resummation predictions for the jet energy profiles with R = 0.7 to Tevatron CDF data in various Pt 
intervals. The NLO predictions denoted by the dotted curves are also displayed. 



We then convolute the light-quark and gluon jet energy functions with the constituent cross sections of the LO 
partonic subprocess and CTEQ6L PDFs 50] at certain collider energy. The predictions are directly compared with 
experiment data, such as the Tevatron CDF data [IB] using the midpoint jet algorithm (53|, as shown in Fig. |9] The 
band represents the theoretical uncertainty caused by the variation of the parameters from C\ = C2 — exp(7s) 1.78 
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to C% = C2 = exp(— 7b) w 0.56, which serves as an estimate of the subleading logarithmic effect that is not included in 
our formula. It is evident that the resummation predictions agree well with the data in all Pt intervals. Although there 
is slight difference between the data and the central values of the resummation predictions, the deviation is within the 
theoretical uncertainty. The NLO predictions derived from (1, Pt, v\, R, r) are also displayed for comparison, 
which obviously overshoot the data. The resummation predictions for the jet energy profiles are compared with 
the LHC CMS data at 7 TeV [13] from the anti-kt jet algorithm [H3| in Fig. QUI which are also consistent with the 
data in various Pt intervals. Since we can separate the contributions from the light-quark jet and the gluon jet, 
the comparison with the CDF and CMS data implies that high-energy (low-energy) jets are mainly composed of 
the light-quark (gluon) jets. It indicates that our resummation formula has captured the dominant dynamics in a jet 
energy profile. Hence, a precise measurement of the jet energy profile as a function of jet transverse momentum can be 
used to experimentally test the production mechanism of jets in association with other particles, such as electroweak 
gauge bosons, top quarks and Higgs bosons. 
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FIG. 10: Resummation predictions for the jet energy profiles with R — 0.7 compared to LHC CMS data in various Pt intervals. 
The NLO predictions denoted by the dotted curves are also displayed. 



A careful look at Figs. H] and [TU] reveals that the resummation predictions fall a bit below the data, as the jet 
transverse momentum Pt increases. One of the reasons for this deviation may be traced back to the kinematic 
constraint for the real soft gluon emitted from the special vertex in Eq. (l44l) . This constraint will include too much 
radiation outside the inner jet cone r into the estimate of the energy profile, especially when the jet axis of the rest of 
particles moves toward the edge of the inner jet cone. The extra radiation can be regarded as a power correction to 
the energy profile in the small r region, because its effect is proportional to r. Since more radiation will be included 
as r increases, the energy profile at large r has been overestimated in our formalism. The energy profile is normalized 
to unity at r — R, so the overestimate actually causes suppression of the distribution at small r, explaining the little 
falloff of the resummation predictions in comparison with the data. When Pt grows, the power correction in the small 
r region is strengthened due to the narrowness of the jet, explaining why the deviation becomes more obvious at high 
Pt- The above reasoning suggests a more restricted phase space for the real soft gluon in order to reduce the power 
correction and to improve the consistency between the predictions and the data. This subject will be investigated 
in a future work. Besides, we note that the effects from hadronization and underlying events on jet energy profiles 
have been estimated by using the PYTHIA code and removed from the published Tevatron CDF data [2g. On the 
contrary, these effects have not been removed in the published LHC CMS data 



V. CONCLUSION 



We have developed a theoretical framework for studying jet physics based on the QCD resummation technique in 
this paper. The evolution equations for a light-quark jet function and for a gluon jet function have been derived and 
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numerically solved in the Mellin (N) space. The inverse Mellin transformation from the N space to the jet mass 
space was performed, which demands the inclusion of the nonperturbative contribution in the large N region, in order 
to avoid the Landau pole, and to phenomenologically parameterize the effects from hadronization and underlying 
events. It has been observed that the nonperturbative contribution is crucial for describing the jet mass distribution 
in the low invariant mass region. The needed nonperturbative parameters were determined by fits of the resummation 
formula including the nonperturbative contribution to the PYTHIA predictions for the light-quark and gluon jet 
distributions at certain jet momentum and cone size, which were then employed to make predictions for other kinematic 
configurations. The above complete resummation formula, convoluted with the LO partonic hard scattering matrix 
elements and PDFs, have led to the jet mass distributions in good agreement with the Tevatron CDF data at different 
jet momenta and cone sizes. Our solutions for the light-particle jet functions are ready to be implemented into 
factorization formulas for jet production cross sections from various processes. 

We have also derived the evolution equations for the light-quark and gluon jet energy functions. With the jet 
invariant mass being integrated out, the evolution equations can be straightforwardly solved in the Mellin space. 
The energy profiles were then predicted by convoluting the solutions with LO partonic hard scattering and PDFs. 
It has been checked that the resummation results for the energy profiles associated with a light-quark jet and a 
gluon jet agree with the PYTHIA simulations. We have demonstrated that the resummation predictions for the jet 
energy profiles are consistent with the Tevatron CDF data and the LHC CMS data within the theoretical uncertainty, 
while the NLO predictions overshoot the data. It should be emphasized that our formula for this jet substructure is 
insensitive to the nonperturbative contribution, and does not involve tunable parameters. Hence, the agreement with 
the data is a highly nontrivial success of the perturbative QCD theory. Besides, an improvement to reduce the power 
corrections to the predicted energy profiles can be done and will be investigated in a forthcoming paper. 

Since final states observed in experiments are usually composed of quark and gluon jets, jet substructures are 
sensitive to the ratios between quark and gluon contributions in a given kinematic region. It is also known that 
the components of the quark and gluon jets are related to the initial-state PDFs. For example, the quark (gluon) 
jet component can be related to the initial-state gluon (quark) PDF in the W boson and jet associated production. 
By analyzing the ratio between the quark and gluon contributions to jet substructures, we may extract additional 
information on the PDFs, especially on the gluon PDF in the small x region. On the other hand, new physics beyond 
the SM introduces more hard subprocesses, which may contribute differently to quark and gluon productions in final 
states. Therefore, a jet substructure, e.g., the jet energy profile, can be used to search indirectly for new physics in 
the region, where PDFs are relatively stable, when both theoretical predictions and experiment data become precise 
enough. 

At last, we reiterate that our framework is ready for the extension to the study of heavy-particle jets produced 
at the LHC, which contain energetic light decay products. For instance, a boosted top quark at the TeV scale will 
appear as an energetic jet, when it decays through its hadronic modes. Likewise, a boosted W, Z, or Higgs boson 
decaying into jet modes at the TeV scale will also appear as an energetic jet. The heavy-particle jet function and 
energy profile can be defined at a high energy scale in a similar way in the factorization theorem as presented in this 
work. The additional ingredient is the factorization of the light final states from the heavy-particle jet at the lower 
heavy-particle mass scale, for which the conventional heavy-quark expansion can be implemented. The solutions for 
the light-particle jet functions and energy profiles established in this work will serve as the inputs of this factorization 
formula for the heavy-particle jet. The above illustrations manifest potential and broad applications of our formalism 
to jet physics. 
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Appendix A: NLO JET FUNCTIONS 

In this Appendix we calculate the NLO light-quark and gluon jet functions by expanding Eq. ^ to 0(a s ), and 
demonstrate the cancellation of infrared divergences between the virtual and real corrections in the Mellin space. 
After regularizing the UV divergence in the MS scheme, the NLO virtual correction to the light-quark and gluon jet 
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functions are given by 
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respectively, where A 2 is an infrared regulator, and the Wilson line direction has been chosen as n — (l,n. x ,0, 0) for 
convenience. The quark-loop contributions to the gluon jet function will be elaborated at the end of this Appendix. 
The explicit expressions for the NLO real corrections to the light-quark and gluon jet functions are written as 
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respectively, where the polar angle of the radiated particle momentum has been constrained to be within the cone 
size R. In the Mj — » limit and without restricting the phase space of the soft radiation, i.e., with R — > it, the large 
logarithms in the above expressions are collected into 
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This isolation of the P-independent soft contributions at NLO has followed the treatment of the evolution kernel from 
the real soft gluon emission in Eq. (|20|). 

Combining the NLO real and virtual corrections to the light-quark jet function in the Mellin space, we arrive at an 
infrared finite expression 
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in which the infrared regulator A 2 has disappeared. Those /V-dependent terms suppressed by 1/ N have been dropped, 
whose effect is expected to be minor. Similarly, the NLO gluon jet function is given, in the Mellin space, by 
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Applying the derivative v 2 d/dv 2 in Eq. (|27| to the above expressions, it is easy to see that the double logarithms 
reduce to single logarithms, which contribute to the kernel G + K in Eq. Since the double logarithms are n'- 
independent, G + K is //-independent, and satisfies the RG equations in Eq. ([23]). Choosing the renormalization scale 
p 2 = C 3 2 R 2 P^/(Nu 2 ), the above NLO jet functions become 
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The choice of fi depends on v 2 in the way that we have fi ~ O(RPt) as v 2 — v 2 n = Ci/(C 2 N) for the initial 
conditions, which then do not contain the large logarithms In N. The NLO initial conditions of the Sudakov evolution 
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are derived from Eqs. (|A9[) and (|A10I) . respectively, with Cf 
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in Eq. ([2]) involve the Wilson links on the light cone along the vector £. Setting v 2 
reproduce the In A terms in these original definitions at NLO, leading to the final conditions 
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It is seen that as the integration variable v 2 in Eq. (|2"T)l varies from v 2 n to 2/|, the scale p? varies from 0(R 2 P^) to 
0(R 2 P^/N). The latter describes the soft and collinear radiations in the jet mass distribution appropriately, because 
they mainly occur at a lower scale. 

The NLO terms in the expansion of the Sudakov exponent contain 
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where C/ = Cp or CU, for S q or S^, respectively. Combining the above expansion with Eqs. (|A11|) and ()A12|) . it is 
straightforward to show that the resummed jet functions in Eqs. (|3~T|) and (f3"8")) indeed agree with the final conditions in 
Eqs. (IA13[) and (IA14I) at NLO, respectively. That is, our resummation formalism is matched to the NLO jet functions 
with /1 2 ~ 0(R 2 P^/N), implying that the single logarithm introduced by our choice of p? has been also summed into 
the Sudakov factor. The all-order summation of this single logarithm corresponds to the RG evolution in p? from 
p 2 = C 2 R 2 P% to p 2 = (C|Cii? 2 P|)/(C 2 A). 
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FIG. 11: Contour for the integration variable y in Eqs. (|31|l and (|33|l . 



At last, we discuss the treatment of the virtual and real quark-loop contributions to the gluon jet function 
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respectively. In the Mj — > and i? — >■ 7r limits, Eq. f|A17l) gives 
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and the infrared finite expression 
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With our choice of /i 2 , the final condition from the quark- loop contributions is written as 

J [ > ^W + "V+« jfinal ~3vri? 2 P 2 s ^ C 2 iV A 3 C 2 J' 1 j 

which has been added into Eq. (|38|) . The absence of the logarithm IniV implies that the quark-loop contribution is 
not important, as verified in the numerical analysis. 

The initial conditions of the jet functions, namely, the prefactors of Sf in Eqs. (fM|) and (|3"5|) are evaluated at the hard 
scale /i ~ O(RPt)- After applying the inverse Mellin transformation to obtain the jet functions jJ LL / NLO (M 2 ., Pr, P), 
which is inserted into Eq. (|41l) to obtain theoretical predictions, the hard scale /i ~ O(RPt) remains. Since J^ V + 
j(i)R,mym organized in the resummation formalism, the regular piece — jW fi,asym has to be added back in 

order to reproduce the complete NLO corrections to the jet functions. This piece is also evaluated at the hard scale 
fx ~ O(RPt) in Eq. P0"j) . Similarly, the regular piece of the quark-loop contribution, jjjJ^q — J^\ R 'q Svm , should be 
included too, which has been combined into Jg^ R — j^ R ' asym on the right-hand side of Eq. P0"|) . 

Appendix B: INVERSE MELLIN TRANSFORMATION 



Because the evolution equations were solved in the Mellin space, we need to perform the inverse Mellin transfor- 
mation to get the solutions in the space of the jet invariant mass. As stated in Sec. Ill, the argument \i 2 of a s (ii 2 ) 
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in the Sudakov integrals should be treated as a complex number in the inverse Mellin transformation. Besides, the 
argument becomes very small (lower than the QCD scale Aqcd) in the large N region, and the running coupling 
constant suffers the Landau pole problem [43, 55]. To avoid the Landau pole, we introduce a critical scale (i c , below 
which the running coupling constant is frozen to a constant value a s {n 2 ). To be precise, the following prescription is 
proposed 
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We have adopted the perturbative expansion of a s in the numerical analysis 
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The variable AT, appearing in the lower bound of y in Eqs. (|3T|) and (|3"3")l , should be also treated as a complex 
number in the inverse Mellin transformation. The contour in the complex y plane is depicted in Fig. 1111 according to 
which an integral over y is handled in the following way, 
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with the variable changes y x = Ci/\N\ + (C 2 - Cx/\N\)t and y 2 = Ci/\N\ exp(iArg(l/AT)(l - t)). 




FIG. 12: Conventional contour of N adopted in inverse Mellin transformation. 




FIG. 13: Contour of iV adopted in our inverse Mellin transformation. 
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The inverse Mellin transformation for the jet function is defined as 
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with C labelling a contour of N. The conventional contour of N shown in Fig. [T^] is not suitable for a numerical 
approach using a grid file, since different jet masses require different parameterizations of this contour in order to get 
enough information in the large N region. Instead, we choose the contour depicted in Fig. 1131 The inverse Mellin 
transformation along the upper-half part of this contour is written as 
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The expression associated with the lower- half contour can be obtained by taking the complex conjugate of Eq. (IB6|) . 
The parameters involved in the integral variables N% and N 2 are set to c = 5, L = 10, and e = 10~ 6 in our numerical 
analysis. 



Appendix C: NLO JET ENERGY PROFILES 

In this Appendix we calculate the NLO light-quark and gluon jet energy functions defined in Eq. (142]) . Due to their 
lengthy expressions, we focus only on the logarithmic terms below. The NLO virtual corrections with the factorization 
scale fi — Pt and the real corrections to the Mellin-transformed jet energy functions are given by 
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respectively. Combining the virtual and real corrections, we derive the infrared finite NLO expressions 
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in the r — > limit, where the infrared regulator A 2 has disappeared. 

The singular NLO terms of the resumed jet energy functions in Eq. (|56p are given by 
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Substituting the vector n fi = (1,(4 - i? 2 )/(4 + i? 2 ),0,0) for n in Eqs. (IC5[) and (IC6[) to obtain the final con- 
ditions, and choosing the 0(1) constant C = exp(5/2) (C = exp(17/6)) for the light-quark (gluon) jet in 
Eq. (|C7|) . we observe the consistency between Eqs. (|C5j) and ()C6|) . and Eq. (|C7|) . That is, the resummation for- 
mula in Eq. ([55]) has collected the important logarithms in the NLO jet energy functions. The complete expres- 
sions for the NLO initial conditions corresponding to the choice n in — (1, (4C| — r 2 C 2 )/(4Cf -I- r 2 C 2 ),0,0), and 
their convolution formulas with the LO partonic subprocesses and PDFs, can be downloaded from the web site 
http://hep.pa.msu.edu/people/yuan/public_codes/JETENPRO/code_energy_convolute_public.zipl 
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